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Lesions in central nervous system (CNS) and their growth leads to debilitating diseases like 
Multiple Sclerosis (MS), Alzheimer's etc. We developed a model earlier [6] [10] which shows how 
the lesion growth can be arrested through a beneficial auto-immune mechanism. The success of the 
approach depends on a set of control parameters and their phase space was shown to have a smooth 
manifold separating the uncontrolled lesion growth region from the controlled. Here we show that 
an optimal set of parameter values exist which minimizes system damage while achieving control of 
lesion growth. 

PACS numbers: PACS numbers: 89.75.-k, 89. 75. He, 87.19.Xx, 87.19.La, 87.18.-h, 05.90.+m 



I. INTRODUCTION 

Human brain normally contains more than eight bil- 
lion neurons and, with their respective axonal connec- 
tions, create a huge complex network. Multiple sclerosis 
(MS) is a disease of active demyelination of neurons with 
destruction of brain functionality [1 . MS directly affects 
the connectivity of the affected region by disrupting sig- 
nal transmission but the pathogenetic mechanism of MS 
is still unclear. Various clinical studies have been at- 
tempting to explain the pathophysiology of the disease. 
Some studies proposed MS as a neurodegenerative disor- 
der where our own immune system, fooled perhaps by 
molecular mimicry, causes destruction of own nervous 
system [IT] . 

Clinical pathology studies have demonstrated hetero- 
geneity of MS in the immunopathological profiles of the 
lesions in different cases. Lesions are categorized as au- 
toimmune encephalomyelitis (patterns I and II) or oligo- 
dendrocyte dystrophy (patterns III and IV) using a broad 
spectrum of immunological and neurological markers [5] . 
In MS with pattern III lesions, demyelination is induced 
by a functional loss of oligodendrocytes, possibly as a re- 
sult of infection with an unknown virus. European Union 
FP6 NeuroproMiSe project study provides the evidence 
for an association between Epstein virus infection and 
MS, but unequivocal proof of the same is still lacking 

In our previous work, we invoked apoptosis initiated in 
peripheral regions of infection by an intercellular distress 
signaling mechanism to block pathogen advance. This 
model resembles the action of firemen who control wild 
fire by burning peripheral vegetation. Our model demon- 
strated that spread of the infection can be contained by 
such a systematic immune response initiated in periph- 
eral regions of infection [5J [TU] . 

An alternative model supports the possibility of tis- 
sue pre-conditioning in the same peripheral layers as a 
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method to control the radial spread of lesions [jj]. In 
such a mechanism, stress signals emanating from the 
pathogenically affected cells, just as in the apoptosis 
model above, induce expression of heat shock proteins 
and hypoxia-inducible factors in the cells of the rim re- 
gion of lesions. Due to the neuroprotective effects of the 
expressed proteins on the tissues, cells are strengthened 
to resist further infection. It is possible that both of the 
above approaches are employed in the real systems [5]. 

Lesion growth affects structural connectivity of the 
neural network and results in losses in functionalities. We 
adopt graph-theoretic techniques in this paper to quan- 
tify the network robustness and loss of functionalities in 
the complex brain networks. We also establish here that, 
while a range of values of control parameters enables con- 
trol of infection, an optimal set exists, which will mini- 
mize damage to the system. 



II. MODEL 

Extraction of complete axonal projection map of all 
neurons of mammalian brain is currently out of reach 
[T3] and, so, axonal connectivity of the CNS is here rep- 
resented by a fixed radius random graph. This consists of 
nodes generated randomly in a unit square with connec- 
tions established between each node and all its neighbors 
within a fixed radius R. Initially, health status of all edges 
(axons) are assigned an arbitrarily assigned maximum 
value (set to 1 here) and these weights become zero when 
they are killed by infection or apoptosis. Pathogenic pro- 
cess is driven by probabilistic events with a pre-assigned 
probability value of damaging edges in each visit to an 
edge. 

Damages to its edges at a particular node triggers an 
alarm signal when the health status of the node falls be- 
low a threshold, r a ;. This stress signal is propagated to 
all connected nodes from where it is propagated further. 
Apoptosis process is initiated at a node when accumu- 
lated alarm signals at the node reaches a threshold, Tbf- 

Apoptosis affects a circular region with a radius pro- 
portional to a parameter, Cbf, and a higher value of Chf 
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implies apoptosis of a larger region. In apoptosis, a cir- 
cular region around the activated node gets reset to zero. 
No additional signals are generated at these nodes to the 
alarm signals generated in the pathological process. 

Fig. [T] gives snapshots of three different possble scenar- 
ios. In Fig. [T^,, when the control level is low, infection 
is seen to consume the entire system. In Fig. [IJd, level 
of control is still insufficient and this only aids to spread 
the damage further because apoptosis also adds to the 
killed edges. In Fig. [IJ;, we finally have the right level of 
control and infection is well arrested. 

In strengthening process, unlike in apoptosis, weigh- 
tage of edges in the circular region controlled by Cbf is 
increased to resist against infection. We assume here 
that, once strengthened, those edges become immune to 
further infection. A detailed view of our model can be 
had from our previous publications [SJ [TU] ■ 



Network Measures 

The status of the system is encoded into an adjacency 
matrix and quantitative study of loss of functionalities 
due to the infection and associated system protective re- 
actions can be explored through parameters that show 
network quality such as these listed below. 

Flow Coefficient: It is a measure of "local cen- 
trality" to estimate the capacity of a node to conduct 
information flow between its neighboring nodes. It is 
calculated as the number of actual paths of length two 
divided by the number of all possible paths of length two 
that traverse a central node [2]. 

Global Efficiency: Path length is defined the mini- 
mum number of edges that must be traversed to go from 
one node to another. Global efficiency is related to aver- 
age inverse shortest path length and is inversely related 
to characteristic path length. It is numerically easier to 
compute than characteristic path length and is useful to 
estimate topological distance between elements of discon- 
nected graphs [7]. 

Edge Density: It is the proportion of connections 
that exists relative to the number of potential connec- 
tions of a network. An edge density of 1, corresponding 
to a percentage of 100%, would mean that all possible 
edges exist [T2"] . 

Total Wiring Length: It is the total number of walks 
through all possible paths that never visit the same path 
again. One recent study evaluated importance of the 
optimal wiring length in brain network projections which 
makes this measure important for evaluating connectivity 
®. 

Snapshot of the dynamical events are taken at each 
time step and stored in the form of adjacency matrices 
which are then analyzed to estimate values of the above 
parameters. Identifying the quality of the complex net- 
works needs a multi-level approach with different indices 
such as those listed above. 



III. RESULTS 

As in our previous study, health status of the edges are 
either alive (=1) or dead (=0) and probability of patho- 
logic damage was set to be pd = 0.33. Number of nodes 
was set to 400 and pathological damage was initiated at 
the center of a unit square. We aimed to determine the 
total damage to the network for different value of the 
three control parameters (r a /, Tbf and Cbf) under both 
strengthening and apoptosis control methods. 

In our previous studies, we indicated that the patho- 
logical process is always controlled if C& / value lies above 
a critical value, C^{r a i, Tbf). We had also postulated 
that damage to the system will increase with subsequent 
increase in Cbf due to the increased aggressiveness of 
apoptosis (see our Fig. 3 in [TO])- We illustrate this aspect 
below and focus here on the infection-controlled cases for 
a series of different parameter values in (t q ;,T{,/ ,Cbf) to 
identify the optimal values in the parametric space which 
aids us to minimize the damage to the system. 

Fig.[2ja) shows that, if we monitor the total number of 
edges damaged by infection after it has been controlled, 
we find it to progressively reduce as Cbf is increased from 
1 to 3 (lesser values when Tbf is smaller). This continues 
till a level of minimum damage is attained, after which 
no further reduction is possible however much we increase 
Cbf - This is what we had anticipated in an earlier work. 
On the contrary, the total damage to the system increases 
for further increase in Cbf. This is seen in Fig.^b). This 
growth in total damage is due to the unnecessary damage 
caused by the apoptosis process. Increase in Cbf leads to 
larger areas being subjected to apoptosis in each apop- 
totic event and this is clearly unnecessary in this case. 
Fig. |2jc) strengthens our argument. It is seen here that 
the time taken to control infection attains its minimum 
for the same values of Cb / where minimum damage due to 
infection occurs as well as the minimum in total damage 



As we can see in the insets to Fig. [2 Cbf value increases 
approximately linearly with Tbf. This happens because 
when Tbf increases, the reduced sensitivity of the alarm 
signal causes apoptosis process to be initiated with a de- 
lay and, it then takes longer to control the infection. A 
similar situation occurs when Tbf is held fixed and T a i 
is allowed to vary. However, in this case, smaller values 
of T a i require larger values of Cbf to achieve control of 
infection (see [TU]). 

In the tissue preconditioning process, a very similar 
scenario results as in the above. However in this case, 
we do not see an increase in total damage for values of 
Cbf > C^ because the total damage is all due to infec- 
tion alone and, when the infection process is controlled, 
there is no more damage because of the absence of the 
apoptosis process. In biological systems, it is possible 
that a control system to bring to a halt, when needed, 
the apoptosis process, as also the tissue preconditioning 
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FIG. 1. Snapshots of the network for three different cases of lesion dynamics, (a) Uncontrolled infection (red) killed most of 
the healthy edges (green) because insufficient apoptosis process (blue) failed to check it (r a ;=0.5, C(,/=0.1 and T(,/=0.5). (b) 
War between apoptosis process and infection is on, but the intensity of apoptosis process is not enough to check the infection. 
On the contrary, it only succeeds in adding to total damage (r a !=0.5, C&/=0.8, and T6/=0.5). (c) Success of infection control 
attained with minimal damage by the efficient action of apoptosis process (r a i = 0.5, C(,/=1.5 and tj,/=0.5). 




FIG. 2. (a) Variation in number of edges damaged by infection with increase in Cbf values for different Tbf values (indicated in 
figure) and a fixed r a ;=0.5. (b) Variation in total damage (due to edges killed by infection as well as apoptosis) of the system as 
Cbf is increased for fixed r a i= 0.5 and different Tbf values (indicated in figure), (c) Variation in time taken to control infection 
as Cbf is increased for fixed r a i= 0.5 and different Tbf values (indicated in figure). The insets to the figures plot the values of 
Cbf at the minimum of the graphs versus the respective Tbf values. 



process, exists which will take care of unnecessary dam- 
ages or wastage of energy. In any case, our basic message 
is that an optimal set of parameter values exist which 
minimizes damage and wasteful energy leakage. 



In Fig. [3j the evolution of some parameters indicating 
network quality are shown, through the processes of in- 
fection and the ensuing reaction from the immune sys- 
tem. The values of the parameters are shown for both 
the apoptosis and tissue preconditioning approaches. It 
is seen that the quality of the network is preserved better 
in the tissue strengthening process, as is to be expected, 
than the apoptosis process. 



IV. CONCLUSION 

In an earlier study [TU] , modeling of MS infection 
process, or equivalently just the lesion growth in CNS, 
and its possible control by the immune system, was car- 
ried out employing a fixed radius network and signaling 
mechanisms which initiate apoptosis in regions surround- 
ing the infected regions. Akin to firemen who burn pe- 
ripheral vegetation to control spread of wild fires, the 
apoptosis process was shown to be capable of terminat- 
ing the lesion growth. Our approach showed that autoim- 
munity can work beneficially and need not be considered 
to be a defect in the system. In this we were following 
earlier work by Matzinger et al. and others [5]. 

In this study, we have showed that an optimal set of 
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values exist for the control parameters (t q ;, Tbf, Cbf) 
which minimizes damage to the system while at once 
controlling the infection. We have not tried to model the 
process of stopping the apoptosis process after successful 
control of infection has been achieved. Neither have we 
attempted to model how the system may arrive at the 
optimal values of the control parameters. These aspects 
can be studied after the biological elements which can 
implement the above model of autoimmunity have been 
identified. Currently, implementation of our beneficial 
autoimmune model to the real human brain connectivity 
network is being undertaken to overcome limitations of 
the fixed radius model as also to bring the model into 
a realistic environment where identification of biological 
elements that effect the autoimmune response can be at- 
tempted. 



FIG. 3. Variation in different indices denoting quality of 
network as the dynamics between infection and its control 
through apoptosis/strengthening process progress towards 
the final state of controlled infection. The graph is typical 
across all values in the ranges of the control parameters and, 
is depicted here for r ai =0.6, C /=1.6 and t o /=0.6. 
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